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^ , Abstract. The conserved Manna model with a planar absorbing boundary is studied 

in various space dimensions. We present a heuristic argument that allows one to 
compute the surface critical exponent in one dimension analytically. Moreover, we 
I discuss the mean field limit that is expected to be valid in d > 4 space dimensions and 

' demonstrate how the corresponding partial differential equations can be solved. 
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1. Introduction 
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■ Absorbing pliase transitions are a particular class of nonequilibrium phase transitions 
(see [1-4] for recent reviews). These transitions continue to attract attention since 
^ ■ they have no equilibrium counterparts and occur even in one-dimensional systems. The 
^ ■ stochastic process of directed percolation (DP) is recognized as a paradigmatic example 
of absorbing phase transitions. According to the well-known conjecture by Janssen 
and Grassberger [5,6] a continuous phase transition should belong to the universality 
class of DP if the corresponding model has a scalar order parameter with a single 
absorbing state, evolves by short-range interactions without quenched disorder and has 
no unconventional symmetries and conservation laws. Non-DP behavior is expected to 
occur if one of these conditions is violated. In particular, particle conservation leads to 
an autonomous universality class of nonequilibrium phase transitions which is usually 
referred to as the Manna universality class. In the recent years, the scaling properties of 
the Manna universality class were investigated in several mostly numerical works (see [4] 
and references therein) . The same type of critical behavior occurs in various models such 
as the Manna sandpile model [7] (which was introduced in the context of self-organized 
criticality), the conserved lattice gas (CLG) [8] and the conserved threshold transfer 
process (CTTP) [9]. In addition, the Maslov-Zhang sandpile [10] and the Mohanty- 
Dhar sandpile [11], which were believed to belong to DP, were shown to be a member 
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of the Manna universality class as well [12-14]. Moreover, the Manna universality class 
is of fundamental interest since it connects the critical behavior of absorbing phase 
transitions with the critical state of self-organized critical (SOC) systems [15]. Actually, 
SOC sandpile models can be considered as driven-dissipative versions of (closed) systems 
exhibiting absorbing phase transitions [15]. An intriguing consideration of the interplay 
between absorbing phase transitions and SOC in case of the CTTP can be found in [4] . 

As in equilibrium statistical mechanics, critical systems out of equilibrium respond 
strongly to the presence of boundaries. For systems in the DP class, the problem of 
boundary effects was investigated thoroughly in a series of papers, where numerical 
simulations, field theoretical approaches [16-19], series expansions techniques [20,21], 
as well as density matrix renormalization group methods [22] have been applied. 
These studies focused on two different types of boundary conditions, namely, so-called 
absorbing walls, where activity is set to zero by removing particles, and active walls, 
where activity is artificially held at a high level. Usually absorbing walls generate a non- 
trivial surface-critical behavior which can be described in terms of an additional surface 
critical exponent f3s- On the other hand, the impact of active walls can be explained in 
terms of the standard critical exponents only. In all cases the bulk critical behavior 
remains unaffected, i.e. the bulk scaling behavior still belongs to the ordinary DP 
universality class. Analogous results were obtained for the so-called parity-conserving 
universality class of absorbing phase transitions [19]. 

Such boundary effects may also be studied within the Manna universahty class. 
Compared to DP, the implementation of an absorbing wall in the Manna class is more 
subtle since the dynamics at the boundary has to preserve the global number of particles. 
Furthermore, the existing knowledge about boundary effects in the Manna class is still 
limited and systematic studies started only recently [23]. In the present paper we extend 
these studies, discussing the scaling behavior near the boundary, computing the surface 
exponent exactly in d = 1, and analyzing the mean-field limit which is expected to 
describe the system above the upper critical dimension dc = 4. 

2. Numerical results 

The Manna model is defined on a lattice with L"^ sites, where each site i is associated 
with an integer number rij representing the local number of particles. Sites with < Nc 
{rii > Nc) grains are considered as inactive (active), where Nc is a stability threshold, 
usually Nc — 2. The Manna model evolves by parallel updates by simultaneously 
redistributing all particles at active sites to randomly chosen nearest neighbor sites, 
which in turn may become active. Note that this update procedure conserves the total 
number of particles N if periodic boundary conditions are applied. Each parallel update 
counts as a single time step. The Manna model is known to display a continuous phase 
transition from a fiuctuating active phase into frozen configurations (lacking of active 
sites), where the particle density = N/L'^ plays the role of a control parameter. The 
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Figure 1. Conserved Manna model in two dimensions with an absorbing boundary 
in form of a line: Influence of an absorbing boundary on the density of active sites in 
a linear (left) and a logarithmic (right) plot. Parameters: Lattice with N — 200 x 200 
and periodic boundary conditions simulated a,t (j) = (f>c and averaged over 65536 runs. 
The corresponding results in 1 and 3 dimensions are qualitatively similar. 



phase transition takes place at the critical threshold 0c which depends on the spatial 
dimension, the lattice structure, as well as on the particular value of Nc [4]. In order 
to preserve the number of particles, the absorbing wall is implemented by removing all 
particles from active boundary sites and adding them to randomly chosen sites in the 
bulk of the system. 

First we studied the spatio-temporal impact of absorbing boundaries in different 
space dimensions for systems starting with a homogeneously active state. Fig. [T] shows 
the decay of the density of active sites as a function of the perpendicular distance to the 
wall at different times. As expected, the absorbing wall creates a depletion zone that 
grows with time as t^^^. The density of active sites inside the depletion zone scales as 



where x denotes the distance to the boundary and Hs is a new exponent. According to 
the standard scaling theory for surface-critical phenomena this suggest the scaling form 



P 



for X <^t 



(1) 




(2) 



where i? is a scaling function with the asymptotic behavior 




for ^ < 1 
for ^ > 1 



(3) 



Our numerical simulations lead to the estimates 




1.00(2) for = 1 
1.10(2) ford = 2 
1.40(15) for = 3 



(4) 
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The asymptotic behavior of R imphes that the activity ps next to the boundary, e.g. at 
a site adjacent to the wall, decays as 

Ps{t) t-'^^ . (5) 

Using the scaling form ([2]) with the asymptotic behavior ([3]) for ^ ^ 1 one is led to 

pM = p(l,t) = t-"{^y = t"^-"-'!' (6) 
which implies the scaling relation 

= a + Kgjz . (7) 
Inserting the known estimates from Ref. [4] this scaling relation yields the estimates 

(0.86(4) for = 1 
1.14(3) for = 2 . (8) 
1.51(9) for d = 3 

In one dimension our result is in very good agreement with the estimate = 0.85(2) 
reported in Ref. [14]. 



3. Heuristic proof for /t^ = 1 in one dimension 

The numerical result = 1.00(2) suggests that Kg might be exactly equal to 1. This 
conjecture can be supported as follows. It is argued in [24] that the critical behavior 
of the Manna model can be described effectively by two mutually coupled Langevin 
equations of the formQ 

dtp{^, t)= - rp{^, t) - bp\^, t) + VV(x, t) 

+ u;p(x, t)0(x, t) + a^J p{y., t) ?7(x, t) (9) 
at0(x,t) = DVV(x,t). (10) 

Here p(x, t) describes the coarse-grained density of active sites while 0(x, t) represents 
the coarse-grained background density of all particles. As usual, ri{-x,t) denotes an 
uncorrelated white Gaussian noise. Note that Eq. iQ differs from the Langevin equation 
for DP by an additional term ujp(f) which couples the two fields whereas Eq. flTUl) describes 
diffusive reordering of the background field during toppling. Furthermore note that in 
a system with periodic boundary conditions the structure of Eq. ( fTOl) ensures global 
conservation of the background field since dt J d'^x(j){'x,t) = D J d'^a; V^p(x, t) = 0. 

X In our opinion these Langevin equations should be considered as an educated guess. Derivations 
from microscopic models are available [14] but they involve various approximations which are hard to 
verify. 
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In the following we show that the structure of these Langevin equations implies 
Kg = 1 in one dimension. The idea is that the removal of active particles at the boundary 
and its redistribution in the (ideally infinite) bulk of the system depletes not only in 
the coarse-grained activity p{x,t) but also the background density t). This is in 
fact the essential feature that makes boundary phenomena in the Manna class different 
from those of DP. Using the Langevin equation we compute the integrated loss of the 
background density and show that the results in one dimension would be inconsistent 
unless = 1. 

To this end we choose a particular site at distance xq and monitor the loss of the 
background density at this point. Since the depletion zone grows as t^^^, this density 
will begin to decrease at about a typical time to which scales as Xq. How much of the 
background density will be lost between t = to and t = oo? According to Eq. ffTOl) . the 
total loss A0 = 0(xo, to) — 0(a^o, oo) is given by the integral 

Act) = - dt^ixo,t)dt ^ - Ddlp{xo,t)dt. (11) 

J to J to 

Inserting the scaling form ([2]) and substituting the scale-invariant ratio 

^ = (12) 



we obtain 



POO 

-D t-2/^-"i?"(xo/t'/") dt 

J to 

- Dz x^^^-^ ! ° ^1--+-° i?"(^) de , (13) 







where ^o = XQ/t\l^ . Note that this expression is positive since R{C) has a negative 
curvature. 

Since the total loss A0 is a finite quantity, the integral ( |T3l) has to be finite as well. 
Obviously, the integral is finite if and only if the integrand diverges slower than 1/^ as 
^ ^ 0. In this limit the scaling function can be approximated by R{C) ~ i'^" so that 
the integrand diverges as 

e^-"+""i?"(0 - k,{k,-1)C^-'^"^+^' . (14) 

For this reason the integral is finite if one of the following three conditions holds: 

• = or 

• Kg = 1 or 

• as = a + Kgjz > 1. 



The first solution = is unphysical because there would be no depletion zone, 
meaning that we are left with the other two possibilities. Which of them applies depends 
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on the dimensionality of the system: In one dimension the exponent = 0.86(4) is 
clearly smaller than 1, hence the above integral is finite if and only if = 1. In higher 
dimensions, however, as is larger than 1 so that Hs can take arbitrary values. This 
explains why = 1 in one-dimensional systems only, refiecting particle conservation. 



4. Mean field approximation 

In d > 4 dimensions, i.e. above the so-called upper critical dimension dc = 4, the Manna 
model at criticality is expected to display mean field behavior on large scales which is 
usually characterized by simple values of the exponents. In the following we compute 
the surface exponent a.^ in the mean field limit, verify the validity of the scaling form ([2]), 
and show how the corresponding partial differential equations can be solved numerically. 

As usual, we assume that the mean-field equations have the same form as the 
Langevin equations (!9l)-( fT0l) if the noise term is neglected. Since the boundary breaks 
translational invariance, space-dependence has to be retained. Moreover, the time scale 
of the process may be fixed by choosing 6 = 1. Furthermore, the parameter D can be 
absorbed in u via rescaling of and thus we may set D = 1. The resulting set of mean 
field equations reads 

9tp(x, t) = - rp{^, t) - p2(x, t) + VV(x, t) + cup(x, t)0(x, t), (15) 

ai0(x,t) = VV(x,t). (16) 

If the boundary has the form of a — 1-dimensional hyperplane it is clear that the mean 
field densities will only depend on the distance perpendicular to the wall. Denoting this 
distance by x we may replace the gradient by an ordinary derivative perpendicular to 
the wall, leading to 

dtp= -p' + dlp+iuj<f)-r)p, (17) 
dt<P = dlp, (18) 

where we omitted the arguments x, t. Obviously the critical background field is 
0c = r/u. We iterated these partial differential equations numerically by means of 
standard methods (see Fig. [2]). Because of the scaling form ([2]) 

p{x,t) ^ r"i?(0 (19) 

with the scale-invariant variable ^ = x/t^^^ one expects a data collapse if t°'p{x,t) is 
plotted against ^. As shown in the left panel of Fig. [2l a reasonable collaps is obtained 
for a = 1 and z = 2, the well-known bulk exponents of DP in the mean field limit. 
Similarly, the background density should obey a scaling form 

0(x,t) ~ 0c-r"'i^(O (20) 

where the numerical integration (see right panel of Fig. [2]) leads to the value a' = 1. This 
is reasonable since the first and the third term in Eq. (fT7|) should be equally relevant, 
hence a = a' in the mean field limit. 
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Figure 2. Numerical integration of the mean field equations ([17]) and ((T8)) at the 
critical point. The left and the right panel show results for the density of active sites 
p{x, t) and the deviation of the background density of particles (j){x, t) from its critical 
values (t)c measured at t = 1, 2, 4, 8, ... , 65536. The data is rescaled according to the 
expected scaling forms while the insets display the raw data. 



Although the data collapses become gradually better for large t, they are not really 
perfect, raising the question whether slightly different critical exponents would lead to 
a better data collapse. However, this would be not consistent with the structure of the 
mean field equations (ITTI) and (JTSi) . To see this we eliminate the parameter r by a shift 
of ip + r/uj, solve the first differential equation for and insert it into the second 
one, arriving at a single partial differential equation for the density of active sites p 



{d,p){dlp - d,p) = p\ujdlp-d,p) + p{dld,p-dlp) 



(21) 



Inserting the scaling form (fT9|) one obtains an autonomous ordinary differential equation 
for R{C) (meaning that it does no longer depend explicitely on x and t) if and only if 
a = 1 and z = 2. This leads to a non-linear third-order differential equation of the form 

1 



2e 



2R{^)(2ujR!\i) + iR!{0-2 



+ e^(2/?"(0 + W0 



(22) 



The trivial solution i? = is unphysical while the solution R = 1 describes the 
homogeneous case without boundary. With an absorbing wall one has to impose a 
Dirichlet boundary condition i?(0) = together with the bulk limit R{oo) = 1. The 
second term ^ '■^^^^^ ^'^^ in the second line of Eq. fl22|) requires that either R'{0) or R"{0) 
vanishes, hence the solution R{^) is either linear or quadratic at the origin. The linear 
solution can be excluded because it would imply that the loss of the background density 
according to the second Langevin equation (ITOl) would be zero. Therefore, the solution 
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behaves as -R(0 oc for small ^, hence we obtain the exponent Kg = 2. 

We did not succeed in solving Eq. fl22|) analytically. However, Eq. fl22l) can be 
iterated numerically by setting R{0) = R'{0) = and tuning _R"(0) ^ 0.283272 in such 
a way that the iteration approaches -R(0 = 1 large ^. The data (not shown here) is 
in perfect agreement with the directly iterated mean field equations. 



5. Conclusions 



In the present paper we have investigated the properties of an absorbing hyperplane in 
a critical Manna model. All results are consistent with earlier findings by Bonachela et 
al [14,23] in Id, who used such boundary effects as a criterion to discriminate between 
directed-percolation and Manna scaling. We have extended these numerical studies to 
higher dimensions and have confirmed that the observed scaling properties are consistent 
with the mean field limit. Moreover we have shown that in one spatial dimension, the 
structure of the Langevin equations, which are believed to describe Manna-type critical 
behavior on a coarse-grained scale, allows one to calculate the surface-critical exponents 
by means of a scaling relation. 
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